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Abstract 

Use of density-functional theory in a Aself-consistent 
field framework result in both the ground- and two low- 
est electronicly excited states of the NaCl and LiCl. 
The accuracy of this method is confirmed using a multi- 
configuration self-consistent field method to obtain the 
same states. The overall good agreement between the 
calculated ground and excited potential-energy surfaces 
speaks promising for the computationally simple Aself- 
consistent field method. 
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I. INTRODUCTION 

The development of powerful computers has allowed 
available ab initio methods to calculate electronic ex- 
citations within considerably larger systems than ever 
before^. However, the need to calculate electronically ex- 
cited states within systems of such a size, that the above 
methods are just not practical, is ever present. 

The less computational expansive density-functional 
theory^ (DFT) as proved its value as a theoretical 
method over the years. Typically this method address 
only ground-state properties. However, over the years 
there has been a number of studies extending DFT into 
the realm of electronic excitations. Time-dependent 
DFTi^, GW"*, perturbation-DFT5 and embedded DFT^ 
methods, all show an impressive and promissing progress 
in this field. 

In this paper, DFT within a Aself-consistent field 
(SCF) framework is used to calculate the ground and 
two lowest electronically excited states of NaCl and LiCl. 
The DFT-based ASCF method, which recently has re- 
ceived some justification^ is an extraordinary simple 
method to calculate electronic excited states. Here, state- 
of-the-art multi-configuration self-consistent field (MC- 
SCF) method^ is used to confirm the accuracy of the 
above method by calculating the same electronic excita- 
tions. The calculated PES's obtained by both methods 
are compared both with other theoretical results'' and 
experimentalifi ones. The overall agreement between the 
MCSCF method and the DFT-based ASCF method is 
promising for use of the ASCF method on systems of 
such size that are beyond the present realm of the MC- 
SCF method.ii 



As an extantion, the quantum dynamics of a simu- 
lated photodissociation process is resolved using timede- 
pendent wavepacket propagation, on the obtained PES's. 
The wavepacket is simultaneously propagated on the cou- 
pled diabatic PES's and the distribution of the ampli- 
tudes are monitored at each timestep so that the disso- 
ciation fraction can be determined and also the interme- 
diate dynamics of the photoreaction. 

The organization of the paper is as follows. In section 
mthe model and the computational methods used in this 
paper are described. Results from the calculation are 
presented in section llTll Conclusions are given in section 
IIVI Finally, the appendix state the results of a simulated 
photodissociation of the two alkali-halide systems 



II. THEORY AND COMPUTATIONAL 
METHOD 

The main advantage with an ab initio calculation, such 
as the MCSCF method, is the predictive power in esti- 
mating the fundamental forces that act on the involved 
nuclei, both in the ground state and in an excited state. 
Unfortunately, the highly accurate ab initio method is 
hard to apply to larger systems, due to the high compu- 
tational cost. On the other hand, first-principle DFT 
has proven to be an essential tool in describing large 
systems, whereas it has been restricted to ground-state 
properties so far. However, recently ordinary DFT has 
been extended to include electronically excited states in 
a ASCF-fashion,^ with a working accuracy.^^ 

The optimized geometry was obtained by a restricted 
Hartree-Fock (RHF) calculation using GAMESS Here 
the basis set uses Slater- type orbitals (STO) together 
with 6 Gaussians, hence ST0-6G. Both dimers have a 
C2v point group symmetry, in which the Hamiltonian 
matrix is constructed in a basis that transforms accord- 
ing to the Ai(E+, A),Bi(H), and ^2(2", A) irreducible 
representations. Here only the S"^ symmetry is consid- 
ered since one of the aims is to simulate a photodissocia- 
tion event. In addition, the use of symmetry reduces the 
complexity of the calculation substantially. The wave- 
function obtained from the RHF calculation is used as an 
input to the MCSCF calculation, where 10 electrons are 
distributed in 14 active orbitals. This generate approx- 
imately one million sets of configuration determinants 
for the Ai symmetry group. The MCSCF determines 
the PES's for the ground and the excited ,^ E+ 
states with high accuracy. 

In analogy with the ASCF^^ method to calculate elec- 
tronically excited states within the Hartree-Fock (HE) 
approximation, the DFT-based ASCF method uses dif- 
ferent electronic configuration in the Kohn-Sham (KS) 
model system to represent these electronically excited 
states. Such an application of the KS formalism has for 
long been without any formal justification, but in a re- 
cent article by A. Gorling^ this method retrieves some 
justification and indeed can be viewed as an approxima- 
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tive method to calculate electronic excitations within the 
considered system. 

The methodology of the used DFT-based ASCF 
method is presented in length elsewhere-'^^ and here only a 
short description is given. The basic ingredients are three 
concepts: (i) interpretation of the KS-orbitals in a molec- 
ular orbital (MO) scheme, (ii) discretization of these or- 
bitals and their energy levels, using supercell calculations 
with periodic boundary conditions, and (iii) introduction 
of electron-hole (e-h) pairs in the system, which is equiv- 
alent to an internal charge transfer in the supercell. First 
an ordinary DFT calculation is performed to obtain the 
ground state PES of the system and the KS orbitals with 
discrete energy levels. Then the relevant KS orbitals for 
the desired internal charge transfer process are identified 
as the ones that should be occupied in the ground state 
but unoccupied in the excited electronic configuration. 
Next a hole is introduced in one of these identified oc- 
cupied KS orbitals together with an extra electron onto 
another one that is introduced into the excited config- 
uration. In this way a KS determinant for the desired 
excitation is constructed. Finally, this KS determinant 
is optimized in a self-consistent-ficld calculation, and its 
energy is evaluated as in a normal DFT calculation. The 
total energy difference between the excited- and ground- 
state electronic configurations is identified as the excita- 
tion energy. 

This method is straightforward, when calculating the 
covalent triplet '^S state of the NaCl and LiCl dimers, 
since it can be constructed using only one KS deter- 
minant. However, in the photodissociation process the 
transitions from the ionic ground-state ^ S to the excited 
state ■^E is forbidden. Therefore we apply the so called 
"sum- method"— to calculate the singlet state of the 
NaCl and LiCl dimers. Simply, it means that a weighted 
sum of determinants is constructed, including both the 
■^S and state, and its energy is found by means of 
ordinary minimization procedure. After that, the energy 
for the singlet state can be extracted since the sole energy 
for the triplet state can be found by similar means. 

The first principle calculations presented in this paper 
are performed by means of the plane wave pseudopoten- 
tial code DACAPC^'^. The generalized-gradient (GGA) ap- 
proximationiSiiLi^ is used for the exchange-correlation 
energy-density functional. The wave functions are ex- 
panded in a plane- wave basis set, and the electron- ion 
interactions are described by ultrasoft pseudopotentials. 
The electronic density between iterations is updated by 
means of a Pulay-mixing algorithm. The occupation 
numbers are updated using a developed technique based 
on minimization of the free energy functional. All calcu- 
lations are performed allowing for spin polarization. The 
dimer cases LiCl and NaCl, are calculated using a super- 
cell of a volume of 20x20x20 A^. The reason for using 
such a big supercell is to minimize any artificial effect 
from the periodicity. 



III. RESULT AND DISCUSSION 

Here the results for total energies and aspects of the 
electronic structure, such as the charge density, are pre- 
sented. The electronic ground states for the NaCl and 
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FIG. 1: The different valence electron density profiles of the 
NaCl molecule through a cut along the molecular axis, (a) 
The ground state with its ionic character, (b) The covalent 
state between Na and CI. 

LiCl dimers at their equilibrium distances are ionic, with 
the Na(Li) and CI atoms being positively and negatively 
charged, respectively. However, at infinite separation the 
ground state is a neutral "covalent" configuration, where 
the ionic configuration lies about 1.4(1.7) eV above the 
covalent PES. This value comes from the difference in 
the ionization energy for the Na(Li) atom and the affin- 
ity energy of the CI atom. As the Na(Li) and CI atom 
move toward each other, there occurs at some intermedi- 
ate separation (~ 10 A)^ approximative the PES's for the 
covalent and ionic states have comparable energies, which 
makes an internal charge transfer of the unpaired 3(2)s 
electron on the Na(Li) to the electronegative CI atom 
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and the electrostatic force affect tfie ionic fragment more 
strongly than in the original weak covalent bond. The 
ionic character of the ground state can be seen in Fig.^ 
where the cut through the dimer axis clearly shows how 
the charge is focused around the chloride atom. 

Analyzing the DOS for the system shows that the elec- 
tron transfer moves the 3s-electron from the sodium atom 
to the empty 3p® orbital on the chloride atom, as ex- 
pected. So in order to obtain the excited triplet co- 
valent state, the hole is introduce in the now filled 2>p^ 
orbital and the electron in the parallel spin-channel of 
the 3s orbital of the Na atom. Figure^ shows the co- 
valent character of the excited state with a charge 
density located around the sodium atom, representing 
the 3s electron, and more importantly, a concentration of 
charge between the atoms. Since we are interested in the 
excited singlet ^E state, we construct a combination of 
the triplet state and the singlet state by calculating the 
energy for the KS-determinant constructed by instead 
placing the electron in the anti-parallel spin-channel of 
the 3s orbital of the Na atom. The energy separation of 
the triplet-singlet state is then extracted from this calcu- 
lation. 

In Fig.|3the PES's for the ionic ^E ground state and 
both covalent ■^E,^ E excited states of the NaCl dimer 
are shown, calculated with the MCSCF method and 
DFT-based ASCF method. In Table 1 the bondlength. 
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FIG. 2: The calculated potential energy curves for the NaCl 
system with the MCSCF (left) and ASCF methods (right). 
The ionic ground state is displayed in a solid line while 
the excited ^E,^ S covalent states uses a dashed line. 

and vertical excitation energies calculated with the 



TABLE I: The calculated bondlength for the ionic ^E state 
and vertical excitation energies to both ^E,^ E states, for the 
NaCl and LiCl dimers. Here results from (a) the MCSCF 
method, (b) the ASCF method, are compared with those for 
(c) the valence-bond (VB) method, and experiment. The unit 
is eV. 

NaCl MCSCF ASCF VB Exp. 

Bondlength 2.35 2.32 2.34 2.36 
^E E 5.27 4.92 - 

^E E 5.49 5.17 5.74 5.26 

LiCl 

Bondlength 2.02 2.06 2.07 2.02 

^E E 5.96 5.64 - 

^E E 6.08 5.75 7.28 - 



present methods are compared with values from other 
calculations^ and experiment for both the NaCl and 
LiCl dimers. The overall agreement between the point 
from the MCSCF and ASCF methods is promising. For 
instance, the discrepancy in the (i) vertical excitation for 
a ^E E transition for both the NaCl and LiCl dimers 
is around 0.3 eV, (ii) singlet-triplet splitting is around 
0.05 eV, and (iii) bondlength is around 0.03 A . 



IV. CONCLUSIONS 

The ground and two lowest excited states of NaCl and 
LiCl dimers are calculated. In the paper there are two 
main theoretical methods used; (i) the ab initio MCSCF 
method with high accuracy and (ii) a first-principle DFT- 
based ASCF method with a working accuracy to calcu- 
late both the ionic "'^E ground state and the two ^E,^ E 
covalent states of the systems. 

The overall agreement between results of the MC- 
SCF and DFT-based ASCF methods is very promising. 
For instance, the discrepancy in vertical excitation for a 
^E E transition for both the NaCl and LiCl dimer is 
around 5% and in bondlength around 2% compared to 
each other. 
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VI. APPENDIX 

Knowing the PES's for the system enables the use of 
wavepacket propagation to resolve the quantum dynam- 
ics of the nuclei. However, the quantum dynamics of 
the photodissociation process on its multitude of diabatic 
potentials is a complex problem. As the system makes 
a sudden transition from its ground state to its excited 
state there will be a new set of forces that work on the 
dimer. The potential energy will be released into the 
intermolecular coordinate which might result in a disso- 
ciation of its fragments. The timeevolution of the system 
is determined by the Schrodinger equation 



(6.1) 



where the Hamiltonian matrix for the NaCl dimer looks 
like 



H 



Hcov 

V21 



Vl2 



(6.2) 



Here the diagonal elements Hcov and Hion represent well- 
defined electronic configurations of the NaCl and LiCl 
dimers, i.e. the ionic ground state and the covalent ex- 
cited state. The off-diagonal element V21 is the non- 
adiabatic coupling between the diabatic states, which en- 
ables the wavepacket to bifurcate among the states. It 
has been foundiS that the strength of the non-adiabatic 
coupling term V21 depends exponentially on the position 
of the curve-crossing point Rcross as 



V2I = Vi2ex.p[j Rcross], 



(6.3) 



where for the halide CI atom the parameters are V12 — 20 
eV, 7 = 1.1638 A-\ 

The solution to the Schodinger equation is calculated 
with a timedependent wavepacket method that is based 
on discrete variables and a finite basis representation 
(DVR-FBR),2a The DVR has the advantage that the am- 
plitude of the wavepacket is well defined, and it leads it- 
self to the simultaneous propagation of the wavepacket 
on different diabatic potentials. The propagation of 
the wavepacket is done with the standard split operator 
technique^ where the potential operator together with 
the effect of the non-adiabatic coupling are evaluated in 
the DVR, while the kinetic operator is calculated in the 
FBR. 

As a consequence of the overall good agreement be- 
tween the PES's obtained by the MCSCF and DFT- 
based ASCF methods, the results from the quantum dy- 
namics turns out to be the same within the simulations 
accuracy. Here the quantum dynamics performed by 
the NaCl dimer during photodissociation is found to be 
highly non-adiabatic, with a negligible population of the 
quasi-bound state, whereas the LiCl dimer has a small 
population of the quasi-bound state. It is concluded that 
the population of the quasi-bound state strongly depends 
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FIG. 3: The timeevulotion of |(*(t)|vE'(t))|^, representing the 
probability distribution for the NaCl dimer, at three different 
times. Here the dashed-dotted curve shows the probability 
distribution on the covalent state, whereas the solid curve is 
the probability distribution of the ionic state. 
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FIG. 4: The timeevulotion of |(*(t)|*(t))|^, representing the 
probability distribution for the LiCl dimer, at three different 
times. Here the dashed-dotted curve shows the probability 
distribution on the covalent state, whereas the solid curve is 
the probability distribution of the ionic state. 



on the difference between the ionization potential and the 
affinity for the two atoms involved in the dimer. Hence, 
if one moves up the rows of the alkali-atoms or down the 
rows of the halogen-atoms the population of the quasi- 
bound state is expected to increase substantially. The 
observation of oscillations in pump-probe experiments on 
Nal dimers'^^ should be consistent with this conclusion. 
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